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Abstract 

We study entanglement entropy for parity-violating (time-reversal breaking) quantum 
field theories on in the presence of a domain wall between two distinct parity-odd 
phases. The domain wall hosts a H-1-dimensional conformal field theory (CFT) with 
non-trivial chiral central charge. Such a CFT possesses gravitational anomalies. It 
has been shown recently that, as a consequence, its intrinsic entanglement entropy is 
sensitive to Lorentz boosts around the entangling surface. Here, we show using various 
methods that the entanglement entropy of the three-dimensional bulk theory is also 
sensitive to such boosts owing to parity-violating effects, and that the bulk response to 
a Lorentz boost precisely cancels the contribution coming from the domain wall CFT. 
We argue that this can naturally be interpreted as entanglement inflow (i.e., inflow of 
entanglement entropy analogous to the familiar Callan-Harvey effect) between the bulk 
and the domain-wall, mediated by the low-lying states in the entanglement spectrum. 
These results can be generally applied to 2-|-l-d topological phases of matter that have 
edge theories with gravitational anomalies, and provide a precise connection between 
the gravitational anomaly of the physical edge theory and the low-lying spectrum of 
the entanglement Hamiltonian. 
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1 Introduction and Preliminaries 


In recent years, entanglement entropy has acquired a central position in the study of quan¬ 
tum held theories. It has emerged as a powerful tool to probe the structure of quantum 
states primarily because: (i) it is sufficiently non-local to capture certain global properties, 
and (ii) it is geometric by dehnition and hence universal in its applicability. As a result, en¬ 
tanglement entropy has provided great insights in a wide class of systems such as relativistic 
held theories [1-3], 2-|-l-d topologically ordered phases of matter [4-13], strongly-coupled 
theories with holographic duals, [14,15] etc. There have also been suggestions that entangle¬ 
ment might play a crucial role in understanding the emergence of geometry in the AdS/CFT 
correspondence [16,17]. 

Entanglement entropy is dehned as follows - consider a density matrix p corresponding to a 
pure state dehned on a Cauchy surface S. Let us partition S into two subregions A and A. 
For local quantum held theories, we expect the Hilbert space "Hs to factorize into the tensor 
product l-ij] = I-La ® 'Ha- If this is the case, we can trace over Ha to obtain the reduced 
density matrix 

PA = Trn^{p) ( 1 ) 

which contains all the relevant information pertaining to the subregion A. Then the entan¬ 
glement entropy between A and A is dehned as the von Neumann entropy of pA 

See[pa] = - Tr^A pa) - ( 2 ) 

In this context, the boundary dA of A is referred to as the entangling surface (or entan¬ 
glement cut). It is also useful to dehne the entanglement Hamiltonian (also known as the 
modular Hamiltonian) He Ai terms of pA as 

PA = e-^A (3) 

Entanglement entropy satishes a number of important properties. One such property which 
will be of relevance to us here, is that entanglement entropy is actually a function of the 
domain of dependence D[A] (see Fig. 1).^ This means that if we were to pick a diherent 
Cauchy surface S' and correspondingly a diherent subregion A' such that D[A\ = D[A'\, 
then the entanglement entropy See[M would be equal to See[^] [18]. Of course, the 
entanglement entropy in the vacuum state of a quantum held theory typically diverges, and 
the above property is really satished by a suitably regulated version of the entanglement 
entropy. 

^The domain of dependence D[A\ of the region A is the set of all points p for which every inextendible 
causal curve through p intersects A. 
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D[A] 



A' 


Figure 1: The domain of dependence D[A] of the subregion A (bold line) is the shaded region. If we 
change A to A' in such a way that D[A] = D[A^], then the entanglement entropy remains unchanged. 
The red dots denote the entangling surface dA. 

In this paper, we are interested in the entanglement entropy in the vacuum state of a parity- 


violating (equivalently time-reversal breaking) quantum field theory on three-dimensional 
Minkowski spacetime where we pick the region A to be a half-space. By a parity- 
violating quantum field theory, we mean a quantum field theory which when coupled to 


background gauge or gravitational fields, contains parity-odd terms in its low-energy effective 
action. The quintessential examples of such terms, which actually suffice for our purposes, 
are given by the Chern-Simons terms 



')+■■• ( 4 ) 


where A is a background gauge field, 5 ^ is a background metric and F the corresponding 
Christoffel connection of g. The set of coefficients (t,k, etc., give a characterization of the 
phase of the theory with respect to parity. In this paper, we will mostly be interested in the 
gravitational Chern-Simons term, and consequently we will restrict our attention only to k. 

An interesting scenario in such theories is as follows: if x labels one of the spatial coordinates, 
then consider a situation where the value of the Chern-Simons coefficient k jumps from H- 
to K+ across the codimension-one surface x = 0. We will refer to this configuration, and in 
particular the surface x = 0, as a parity domain wall. Such a domain wall typically hosts 
a 1-1-1-dimensional conformal field theory (CFT), with a non-trivial chiral central charge 
Cl — cr. However, it is well known that CFTs with cl — cr ^ 0 have gravitational anomalies 
- this means that when such a CFT is coupled to a background metric, its stress tensor is 
not conserved [19,20] 



( 5 ) 
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and thus (background) diffeomorphism invariance^ is broken. Here is the consistent 

stress tensor which is obtained by directly varying the CFT partition function with respect to 
the metric, and its non-conservation is called the consistent anomaly. For the parity domain 
wall, the anomaly is remedied by the fact that the bulk theory is also not diffeomorphism 
invariant (in the presence of a domain wall). Indeed, the bulk gravitational Chern-Simons 
term transforms under diffeomorphisms to precisely compensate for the CFT violation of 
diffeomorphism invariance, provided that 




Cl — cr 
487r 


( 6 ) 


In this case, the full theory, i.e., the bulk theory plus domain wall CFT, respects background 
diffeomorphism symmetry. This cancellation of the CFT gravitational anomaly by the bulk 
is an example of Callan-Harvey anomaly inflow [21]. 


Although the above discussion has been somewhat formal, a simple yet concrete example to 
keep in mind is a Dirac fermion on with a mass domain wall, i.e., the mass of the fermion 
is a function m{x) = mQ(p{x) of one of the spatial coordinates x, where ttiq is a positive 
constant and (p{x) smoothly interpolates between —1 and -|-1 as x goes from —oo to -f-cx). It 
is easy to show by solving the Dirac equation that the domain wall traps a left-moving chiral 
(Weyl) fermion in the limit tuq —)■ oo, which can be represented by a free CFT on the domain 
wall with Cl — cr = 1. On the other hand, far from the domain wall on either side, one can 
integrate out the bulk fermion (which recall becomes very massive as tuq —?• oo) to obtain 
the bulk effective action. Indeed, one hnds that this effective action has a gravitational 
Chern-Simons term, where the coefficient precisely satishes equation (6) [22,23]. In fact, 
this is essentially the macroscopic physics of a condensed matter system dubbed a Chern 
insulator [24] which exhibits an integer quantum hall effect (without an external magnetic 
held) due to the time-reversal/parity-violating Dirac mass. In that context, the gapped bulk 
region x < 0 is taken to be the interior of the non-trivial topological Chern insulator, and 
X > 0 is the trivial insulating vacuum. The chiral fermion states on the domain wall are the 
conducting edge states. Another class of condensed matter systems that exhibit the above 
type of anomaly inhow physics are fractional quantum hall states - here the low-energy 
effective held theory describing the bulk is a topological held theory (usually Chern-Simons 
theory), while the domain wall CFT is generically a WZW model [25]. In these systems, 
the gravitational Chern-Simons term in the bulk ehective action is generated by the framing 
anomaly [26]. 


Let us now consider the entanglement entropy for the class of systems discussed above. Let 

^By background diffeomorphism invariance, we mean that the partition function Z[gfj_,^] regarded as a 
functional of the background metric, satisfies the Ward identity + 2d(^^i,)]. 
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= {t, y, x) be coordinates on where t parametrizes time. We take the domain wall to 
lie at a; = 0. We subdivide the spatial slice S at t = 0 into two half-spaces A = {x E T,\y > 
0}, and A = {x E T,\y < 0}, with the surface dA ai y = 0 being the entangling surface (see 
Fig. 2). 



Figure 2: The spatial geometry at t = 0: the darker region x < 0 is a particular parity-odd phase with 
K = K-, while the lighter region x > 0 is a different phase with k = k+. The domain wall x = 0 hosts 
chiral edge states, denoted by bold arrows. The red line y = 0 is the entangling surface, with A being 
the half space y > 0. 

In principle one could consider more general entangling surfaces, but we will focus our 
attention on this case because it simplihes the analysis signihcantly, and will allow us to 
make fairly general statements. Now, we are interested in studying how the entanglement 
entropy between A and A depends on an inhnitesimal Lorentz boost localized around the 
entangling surface with the rapidity rj (see Fig. 3). 

Based on our previous discussion, one expects that the entanglement entropy, being a func¬ 
tion of the domain of dependence of A^ remains invariant under such a boost (which preserves 
D[A\). However, there is a problem here because the entangling surface cuts through the 
domain wall, and the domain wall CFT has a gravitational anomaly. Indeed, it was pointed 
out in [27-32] (and anticipated previously in [33]) that the entanglement entropy of a CFT 
with non-trivial chiral central charge is sensitive to such a local Lorentz boost, and changes 


5 





Figure 3: The boosted Cauchy slice (blue) in the {t,y) plane; the x direction is implicit. Since the 
domain of dependence of ^ (shaded region) remains unchanged, the entanglement entropy is expected 
to remain unchanged. 


by the amount^ 


5See,cft — —f] 


Cl — cr 
24 


( 7 ) 


Thus, for our expectation of invariance of entanglement entropy under transformations which 
preserve D[A] to be true, the entanglement entropy of the bulk must also change under the 
local boost by an amount which precisely cancels the CFT contribution. The primary goal of 
the present paper is to show that this is indeed the case. We will establish this in a number 
of ways. In Section 2.1, we will derive the response of the bulk entanglement entropy to local 
boosts by using perturbative methods developed in [34,35]. In Section 2.2, we will derive the 
same result using the replica trick because this will give us an understanding of the higher 
Renyi entropies as well. Then in Section 2.3, we study the system by conformally mapping 
the Rindler wedge D[A] to M x (where is the two-dimensional hyperbolic space). This 
will give us a way to understand the cancellation between the bulk and the domain wall 
as an inflow of entanglement entropy mediated by the low-lying states in the entanglement 
spectrum. Finally, in Section 2.4 we will present numerical data which demonstrates the 
above physics in the specihc case of a lattice-regularized Chern-insulator model. 


^References [27, 28] use a different set-up involving boosted intervals (instead of boosting locally around 
the entangling surface), and find SSee = —7? ^^ . References [29,30] find 6See for a semi-infinite interval 
to be twice as large as (7), namely 5See = , which might be akin to a covariant form of the 

anomaly. Our result agrees with that of [31,32]. See [30,32] for a more detailed comparison. 
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2 The bulk entanglement entropy 


Let us now delve into the calculation of the change in the bulk entanglement entropy in 
response to a local Lorentz boost around the entangling surface. More precisely, we mean 
the change in the entanglement entropy coming purely from the three dimensional parity- 
violating bulk theory. We will see that the gravitational Chern-Simons term plays a crucial 
role in this calculation. 


2.1 Perturbative method 


In this section, we will use the perturbative techniques developed in [34-37]. Recall that the 
region A we’re interested in is the half-space y > 0, and so the entangling surface is given 
by the line y = 0 (see Fig. 2). We start by performing a local boost in the ^/-direction (see 
Fig. 3) given by the following diffeomorphism 

^(x) = y4>(x) (ydt+ tdy) (8) 

where x'^ = (t, y, x) are coordinates on the full spacetime Here ?7 is a constant equal to 
the local rapidity (or the hyperbolic boost angle) close to the entangling surface, and (j){x) 
is a smooth function such that (j){x) = 1 inside the ball of radius R around x = 0 (with 
R much larger than the thickness of the domain wall, or any short distance cutoffs in the 
theory), and (j){x) —)■ 0 smoothly outside Note that (j) approaches zero at inhnity in all 
directions - this will simplify our calculations by letting us drop boundary terms at inhnity. 
Before moving on, it is useful to Wick rotate to Euclidean signature by dehning t = ite- In 
this case, the above diffeomorphism becomes 


where rjE = —iV- 


= nE4>(x.) (ydr^ - Tedy) 


(9) 


Now, the above diffeomorphism preserves the entangling surface, but changes the metric 
around it 

( 10 ) 

where the Greek indices /i, u now run over the bulk Euclidean indices. The change in entan¬ 
glement entropy to hrst order in the above perturbation is easy to compute and one hnds 
(see Appendix A for more details) 


^iSEEfiulk 


— lim / 


d^xd^,mx) {f^''{x)HE) 

\ / conn, 


( 11 ) 
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where He is the entanglement Hamiltonian, which can be written in terms of the following 
integral over the region H at = 0 

/ oo poo 

dx' / dj/'y'T°°(0, ?/', x') + constant. (12) 

OO J 0 


The correlation function appearing above is the connected Euclidean correlation function on 
and so we may drop the constant term in (12). To regulate this calculation, we have 
cut out a tubular neighborhood of radius a around the entangling surface by restricting the 
integration over x in equation (11) to the region > a? (Fig. 4); we will send a —)■ 0 

at the end of the calculation. Integrating by parts in a;, we hnd 




EE,bulk 


= lim 27ia 


(f d^x f dx' [ dy'y' i'^{x)n''{x) /f^^(a;)foo(0, y', x')) 

^0 /r|+y2=a2 J-oo Jo ^ / 


(13) 


where ni^{x) is the outward pointing unit-normal to the cylinder Te + y"^ = a?- We have 
dropped the term proportional to the divergence of the stress-tensor using the Ward identity 




dy5^{x - x') (fxa{x)'j -F 2a(A - x') (f„)^{x)'j'^ 

d^6{x-x')(TUx)^6xa, (14) 


and the fact that the one-point function of the stress tensor has no parity-odd contributions 
in flat space. We have also dropped all other boundary terms because ^{x) —)■ 0 far away 
from X = 0. Therefore, equation (13) suggests that the change in entanglement entropy can 
be computed from the two-point function of stress tensors. 


Indeed, in theories where parity/time-reversal is violated, there is a parity-odd piece in this 
two-point function [38,39] 

(ff,^{x)fxp{x')) = -2iK e(^(Aa {du)dp) - d^)p)d^) d'^S^ix - x') (15) 

where the coefficient k, is precisely the coefficient of the gravitational Chern-Simons term 
defined in (d).'^ However, recall that in the presence of a domain wall at a; = 0, k is not a 
constant but a function k{x) which jumps across the domain wall, i.e., k{x >> |) = and 
k{x « — |) = K- (where i is the thickness of the domain wall). The diffeomorphism invari¬ 
ance of the full theory (bulk plus domain wall) imposes the constraint (6). Taking this into 

^There is also a term lower order in derivatives of the form [40-43] 

(fp^{x)fxp{x')) = iC e{p.(xadt.)p)d'"5^{x - x') A - 

where C is related to the Hall viscosity response [41,44,45]. It can be checked explicitly that this term does 
not contribute to the integrals in equation (13). 



te 


Z' 

Br 


-5K- 
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Figure 4: A tubular neighborhood of the entangling surface (red dot) is cut out to regulate the 
calculation (the x direction is implicitly out of the page). The dotted line is the cut-off surface 
a?. The entanglement Hamiltonian has support over the blue line. For contact terms in the 2-point 
function, the integral in equation (13) has support only at the point marked with a cross. The shaded 
region is the ball Bji where ^ is supported. 


account, we find that the two-point function of stress tensors has additional contributions^ 


{f^y{x)fxp{x')^ ^ = (^-2iK{x)d'^ + id'^K{x)'^e(p(Xa{du)dp)-6^)pjd‘^)S^{x-x')^ - (16) 

where ■ ■ ■ represents further terms which are not important for the present calculation. Now 
we simply substitute this expression into the formula (13) for and compute the 

integrals. Using the explicit form of ^ given in (9), we hnd 


S^SEE,bulk — 



{.Cl - Cr) 

V —— 


(17) 


where in the second equality we have reverted back to Lorentzian signature, and used equa¬ 
tion (6). Comparing with equation (7), we see that the change in the bulk entanglement 
entropy upon performing a local boost precisely cancels the contribution coming from the 
domain wall CFT. 


2.2 Replica trick 

To demonstrate the robustness of the above result, let us now repeat this calculation using 
the replica trick. This also has the added advantage that it will give us an understanding of 

®These can be obtained by promoting k to a function k{x) and varying the gravitational Chern-Simons 
action twice with respect to the metric. 
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what is happening with higher Renyi entropies. Recall that the nth Renyi entropy is defined 
by 

Sn = ^^HaPa- ( 18 ) 

Note that taking the limit n —)■ 1, approaches the conventional von Neumann entangle¬ 
ment entropy. For integer n, there is an explicit prescription to compute these entropies 
using path-integral techniques [1,15]. One introduces the replica manifold Mn, which is an 
n-sheeted covering of the original (Euclidean) spacetime manifold branched along the region 
A (see Fig. 5). 



Figure 5: The replica manifold for n = 3 has three sheets (the x direction is implicit) which are glued 
along A as indicated - the blue line on sheet one is glued with the blue line on sheet 2 and so on. In 
this case G [0, Ovr) with y? = Gtt identified with (p = 0. As p goes from 27r(/c — 1) to 27rk, we cover 
the kth sheet. 

Then, we have 

Sn,bulk = T~~ Zbuik{Mn) - u In ZbulkiMi)] (19) 

where Zbuik{Mn) is the partition function of the bulk theory on the replica manifold M„, and 
Ml is the original spacetime 

Now we want to ask the same question - what happens to the nth Renyi entropy under the 
local diffeomorphism (9) around the entangling surface? From the above expression, we find 
that the Renyi entropy changes as 

S^Sn,buik = Zbuik{Mn) - n6^\n Z{Mi)j. ( 20 ) 
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We now use the fact that in parity-violating theories, the effective action has a gravitational 
Chern-Simons term 


- In z,uik{Mn) = f ^ a A a t\\ +■■■ ( 21 ) 

where the ellipsis indicates other terms which might be present, but do not contribute to 
the present computation.® Under a coordinate transformation, the Christoffel connection F 
transforms as 

= ( 22 ) 

and consequently the Chern-Simons action transforms as 


5^1n = I [ dKA{dfsC)dT^a (23) 

= / dnACdidf^rf^a)-',: [ d (d^n) A CdT^^^ (24) 

where we have consistently dropped all boundary terms, because ^ 0 at inhnity. The 

integration by parts in going from (23) to (24) is meaningful as long as the geometry of the 
replica manifold is appropriately regulated (see below). The hrst term in (24), in comparison 
with equation (5), makes it clear that the entanglement response to boosts is directly related 
to the consistent gravitational anomaly of the domain wall CFT. The second term in (24) 
vanishes simply because the index (3 has to be along the domain wall for the Christoffel 
symbol to be non-trivial, but then dj^K = 0. Hence, we conclude that 


5gln Znfiulk — —5^hl ZnfiFT 


(25) 


which directly implies that the anomalous response of the bulk entanglement entropy cancels 
the anomalous response of the boundary. Since the anomaly above is the consistent anomaly, 
the gravitational Chern-Simons term - which is not diffeomorphism invariant - is the only 
relevant term in the bulk (see [20,22,42] for a discussion of consistent v. covariant anomalies, 
and their relation to Chern-Simons terms). 


Now all that remains is to substitute the above expression into equation (20). Let us intro¬ 
duce polar coordinates in the {te^v) plane 

r = + (P = tan“^^. (26) 

®For instance, one might worry about the parity-odd Hall viscosity term [42] which was mentioned in 
footnote 4. We now see that this term, though present, cannot contribute to the present computation because 
it is diffeomorphism invariant. 
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With a slight abuse of notation, we take ip G [0,27rn) - when ip crosses 2'Kk (for fc G Z, 
1 < fc < n — 1) we move from the /c-th replica sheet to the {k + l)-th sheet. Now, the full 
n-sheeted replica manifold can be regarded as the space with the metric 

Qn = + n^r'^dd'^ + dx^ (27) 

where we have dehned ip = n9, with 6 G [0, 27r). This metric has a conical singularity along 
the line r = 0 which needs to be regulated. We use the regularized metric [46] (see Fig. 6 
below) 

dna = ^ ? dr"^ + n^r^dO^ + dx^. (28) 

’ 

Clearly, in the limit a —)■ 0 we recover the original cone metric (27). We hnd that the 



Figure 6 : (Left) The geometry of the replica manifold in the {te^v) directions is that of a cone 
embedded in (Right) We can regularize the conical singularity at origin in the (T^,y) plane by 

replacing it with a smoothed-out cone. 


non-trivial Christoffel symbols for the regularized metric are given by 


^ r(r^ + pr ^ aV(n^ - 1) . 

r2 + n2a2 ’ (r2 + a2)(r2 + n2a2) ’ 

After substituting these into equation (24), and using ^ one hnds 

forward calculation 

(1 — 

d^ln Zbulk{Mn) = i7TT]E{K+ - K_)- -. (30) 

n 

Putting this into equation (20) and using ( 6 ), we hnd^ 


(29) 

r 

from a straight- 


X Q - ■ ^L- Cr 

^n^bulk 



(31) 


Finally, taking the limit n —)• 1 and reverting to Lorentzian signature, we obtain 


S^SEE,bulk — V 


{cl - Cr) 

24 


(32) 


^The Renyi entropies in topological field theories are typically n-independent. Here, the n-dependence 
for the bulk entanglement response clearly comes from the presence of the domain wall defect. Indeed, it 
precisely cancels the contribution coming from the domain wall CFT. 
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which indeed matches our calculation in the previous section. Additionally, we repeated the 
above calculation using Cartesian coordinates (instead of polar coordinates) for the replica 
manifold, and we found precisely the same result stated above.® In a way, this calculation 
makes transparent the fact that the cancellation of the anomalous entanglement response of 
the bulk and the domain wall CFT is forced upon us by the diffeomorphism invariance of 
the full theory on replica manifolds. 


2.3 Entanglement inflow 


The calculations in the previous section, while concrete, perhaps hide the physics behind 
the matching of the bulk and domain wall CFT results. We will now describe the bulk¬ 
boundary cancellation of the anomalous entanglement entropy using a more intuitive, albeit 
less rigorous picture of flow of entanglement entropy (or equivalently entanglement ‘energy’) 
which emerges in the Rindler wedge. The Rindler wedge is the domain of dependence of the 
region A (where recall that A in our calculation is the half space ?/ > 0 on the spatial slice 
at f = 0) over which the reduced density matrix is defined (see Fig. 7a). 




Figure 7: (a) A constant x slice of the Rindler wedge. The modular Hamiltonian generates translations 
along X- (b) In Euclidean signature, the modular Hamiltonian generates rotations around the entangling 
surface, i.e., translations in ip. (c) A constant x slice of 5^ x with ^ = ln(z). The r direction is 
compact. The red dashed line represents the cutoff at z = a. The bold blue line represents the region 
A over which the reduced density matrix is defined. 

Therefore, all of the physics relevant to the region A can be equivalently described by going 
to Rindler coordinates 

X = tanh"^ , p = (33) 

®It has been pointed out in [30] that it might be important to treat the coordinate singularities at the 
entangling surface as boundaries, possibly leading to extra contributions. 


13 






















where the transverse coordinate x remains intact. The Minkowski metric in these coordinates 
is given by 

gmndier = -p^dx^ + dp^ + dx^ . (34) 


Importantly, the modular/entanglement Hamiltonian in this pictnre is proportional to the 
generator of boosts, or y translations. It is more convenient to consider the Enclidean version 
of this statement. Consider Wick-rotating from Minkowski space to Enclidean space 
with the metric 

9Eud. = r^dif^ dx^. (35) 


with ip G [0, 27r) being the polar angle in the (r, 93 ) plane (See Fig. 7b). In this case, the 
entanglement Hamiltonian is exactly 27r times the generator of rotations in the (r, p) plane, 
or in other words (^-translations (see Appendix A for a derivation). It is fnrther convenient 
to transition from to the conformally eqnivalent space x (where stands for two- 
dimensional hyperbolic space). Let us use the coordinates® (r, z, x') to parametrize x 
with r G [0, 27r), 2 ; G (0, cx)) and x' G M (Fig. 7c). In these coordinates, the metric on x 
is given by 


— dr"^ + 


dz^ + dx"^ 


(36) 


The map 1 —>■ x in terms of the coordinates {p, r, x) on and the coordinates 

(r, z, x') on x is simply 


T{r,p,x) = p, z{r,p,x)=r, x'{r,p,x)=x. 


(37) 


Since the x coordinate maps trivially to x', we will drop the prime from now on for sim¬ 
plicity. It is clear that (37) is a conformal transformation. This fact has been used quite 
successfully in studying the entanglement entropy of conformal held theories [37,47] . In our 
case, the domain wall theory is indeed a conformal held theory (in the limit of inhnitesimal 
wall thickness), but in the bulk we do not have a GET. However, the object of primary focus 
in the bulk is the parity-violating, gravitational Chern-Simons term in the bulk ehective 
action - this term is indeed conformally invariant (up to boundary terms). Further, in many 
condensed matter systems of interest, such as topological phases of matter, (e.g., topological 
insulators [24,48] and fractional quantum Hall systems [25], etc.) the bulk is ehectively de¬ 
scribed by a topological held theory at low energies, which is certainly conformally invariant, 
while the corresponding edge theories, at hrst approximation, are gapless CFTs. Hence, we 
expect the description in this section to apply to a wide class of such topological systems. 


The reason it is convenient to transition to x is as follows: the region A over which 
the reduced density matrix is dehned maps to a constant-r slice in x The modu- 

®The coordinate r should not be confused with the Euclidean time coordinate te- 
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lar/entanglement Hamiltonian He now generates ‘time translations’ in the r direction, and 
thus can be regarded as 2tt times the physical Hamiltonian on hyperbolic space: 

- /■“ dx dz - 

He = 2ti / — / —Tt-t( 0, a:) + constant. (38) 

J—ao ^ Jo ^ 

Therefore, the reduced density matrix is precisely the thermal density matrix 

with inverse temperature (3 = 27t, and so in hyperbolic space-time, the entanglement entropy 
is nothing but the thermal entropy. Hence, the hrst law of entanglement entropy [17,49] 

6See = SE = 6 (^He) = 5 . (39) 

naturally takes the form of the hrst law of thermodynamics. Therefore, from this point of 
view the anomalous change in the entanglement entropy due to a local boost is equivalent to 
a change in the total energy of the system. This is quite remarkable since it is an explicit map 
between the variation of entanglement entropy and variation of “energy.” Now, due to the 
gravitational anomaly of the domain wall CFT, its intrinsic energy is not conserved, which 
(by the hrst law) is rehected in the anomalous variation of its intrinsic entanglement entropy. 
However, in the full bulk plus domain wall theory, the total energy must be conserved, and 
as such, the bulk must somehow account for the energy lost by the domain wall. We claim 
that this happens as follows: in addition to the chiral states on the domain wall, there are 
also chiral states localized on the (regularized) bulk entangling surface. Being chiral, their 
energy is also not conserved during a Lorentz boost, but in precisely the right way to account 
for the energy non-conservation of the domain wall. Hence there must be anomaly matching 
between the physical domain wall/edge and the regularized entangling surface. 

To see this in more detail, let us hrst focus on the CFT living on the domain wall. The 
domain wall corresponds to the slice of x at x = 0 , which is actually a cylinder with 
the metric 

dz‘^ 

Qdw = dr'^ H- 5 - = dr'^ + (40) 

where in the second equality we have dehned ^ = ln( 2 :). The local boost of the Cauchy 
surface now corresponds to a local translation in the r direction (see Fig. 8 ) 

^ = r]E(j){z,x) dr (41) 

where recall that 0 is a bump function supported over the ball Be around z = x = 0 
(with R much larger than the thickness of the domain wall, or any short distance cut-offs 
in the theory). Ordinarily, this boost would conserve energy (because the Hamiltonian is 
a conserved charge) provided there is no flux through the boundary of hyperbolic space at 
z = 0, which is really inhnitely far away. However, in practice we must regulate this inhnity 


15 


by putting a cutoff at some finite value z = a}^ Then, owing to the chirality of the CFT, 
the flux condition is clearly violated because of the energy-flux through the window between 
r = 0 and t = rj at z = a. Of course, if the theory is non-chiral such that cl = cr then the 
outgoing flux cancels the incoming flux, and the energy remains unchanged. 



e 


Figure 8: The local boost corresponds to a local r-translation in x BI^. The blue line is the original 
constant r surface over which we measure the energy, while the purple line is the new surface after the 
boost. The arrows indicate that in a chiral CFT, we can have an energy-flux through the cutoff surface 
(dashed red line) during such a deformation. 

It is straightforward to quantify how much energy is lost/gained during this deformation. 
We know that the expectation values of the holomorphic and anti-holomorphic stress tensors 
on a cylinder are given by 

{T) cylinder = {T) cylinder = “^ ( 42 ) 

We can then use the above results to compute the flux through infinity, and one finds 

2tt5Edw = y (^Tc('r,ln(a))) = 27ri J dr (^{f) - (T)j = (43) 

where the subscript DW refers to domain wall. We see that the above result is essentially 
proportional to the Casimir momentum density along the domain wall. The hrst law of 
entanglement entropy then gives 

5See,dw - -^Ve —^— - -V —^— (44) 

So we conclude that the change in the domain wall entanglement entropy is caused by energy 
flowing out from the boundary at z = a. 

with all aspects of entanglement calculations, it is essential to regulate the theory. The cutoff at 
z = a maps back to the exclusion of a small cylindrical region around the entangling surface, which is a UV 
cutoff. 
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Let us now revisit the situation from the bulk point of view (see Fig. 9). As we discussed 
above, since UV divergences arise from the inhnite amount of space close to the bulk entan¬ 
gling surface 2 ; = 0, it is necessary to regulate the geometry in this region. Here we choose 
to put a hard cutoff at z = a (which can be thought of as a stretched Rindler horizon), 
analogous to the brick-wall model of black hole thermodynamics [50]; another possibility is 
to use a lattice regularization. In such a regularization, it is expected that there must also 
exist chiral states localized on the cutoff surface z = a, in addition to the chiral states on 
the domain wall x = 0 [8,51-59]. The existence of these states is forced upon us by the re¬ 
quirement of anomaly cancellation for the theory on x (or equivalently on the Rindler 
wedge). In particular, the modes supported in the region x < | of the entangling surface 
must have net chirality cl — cr = 4:8%K-, while those in the region x > | must have chirality 
Cl — cr = 48%K+, where i is the thickness of the domain wall. It is possible to demonstrate 
the existence of these states explicitly for free fermions by solving the Dirac equation on 
M X with brick-wall boundary conditions; this is done in Appendix B. However, we re¬ 
iterate that they are guaranteed by diffeomorphism invariance (i.e., anomaly cancellation). 
These states present a natural resolution to the apparent loss of energy from the domain 
wall - the energy which is lost in the domain wall CFT at z = a is simply transferred to the 
chiral states localized on the regularized entangling surface in the bulk. 

In order to conhrm this, we may approximately describe the chiral edge states localized on 
the regularized entangling surface as a CFT on the cutoff surface z = a. We now compute 
the energy flux through a constant x slice 

dr <^f™(r,x)^ (45) 

along the cylinder at = a, for x < —£/2 and x > Repeating the arguments given 

above involving the Casimir momentum density, we hnd that the change in the energy of 
the region x < — | of the cutoff surface is given by 

i 

6Ees{x < “ 2 ^ ^ -2%ir]EK.- (46) 

while for the region x > we have 

i 

5Ees{x > 2 ) = +‘^'xir]Ei^+ (47) 

where we have used the subscript ES to indicate that these quantities are computed along 
the cutoff entangling surface. Therefore, the total change in the energy during the boost, 

the limit a —>■ 0, the metric on the cutoff surface z = a is conformally equivalent to the degenerate 
metric lim^_,.o . However, the computation of the flux of Casimir momentum is insensitive to 

this degeneracy. 
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Figure 9: The boost from the bulk point of view for the case = 0, / 0: the blue plane is 

the original constant r surface over which we measure the energy, while the purple surface is the new 
surface after the boost. The black and brown arrows indicate the energy-flux associated with edge states 
localized on the domain wall and bulk entangling surface respectively. The energy flux which leaves the 
domain wall enters the bulk along the entangling surface. The direction the edge states flow along the 
entangling surface depends not only on the nature of the domain wall (i.e. (k+ — k-)) but on the way 
the bulk theories are regularized on opposite sides of the domain wall (i.e. the specific values and 
K_). For this figure we have chosen a regularization where k_|_ = 0 and k_ / 0. 

coming from the entangling surface is 

6Ees = 2TTir]E{K+ - /t_) = 

which precisely accounts for the energy-flux out of the domain wall (43). By the hrst law 
of entanglement entropy, we thus conclude that SSee,es = —SSee,dw- Note that this 
result does not depend on much of the detailed properties of the spectrum of modes on the 
entangling surface. 

The chiral states which mediate the inflow of energy flux between the bulk and domain wall 
essentially constitute the low-energy spectrum of the entanglement Hamiltonian He- We 
have noted above that since this is the Hamiltonian in the above Rindler space construction, 
the diffeomorphism invariance of the full theory forces the existence of these chiral states in 
its spectrum. Indeed, in materials exhibiting topological phases of matter, it is known that 
there is a correspondence between the low-lying spectrum of the entanglement Hamiltonian 
and the low-energy spectrum of the physical Hamiltonian in the presence of a physical 
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boundary [8,51-59]. We have seen above that the cancellation of the bulk and boundary 
responses to a local boost can be explained as an inflow of (entanglement) energy associated 
exactly with these low-lying states in the entanglement spectrum, which we interpret as an 
inflow of entanglement entropy. 



Figure 10: Two domain walls with opposite chiralities at a distance d. The arrows indicate energy-flux 
along the domain walls and the bulk entangling surface. In the limit d —)• 0, the walls can be gapped 
out, but we are still left with chiral states localized on the cutoff entangling surface. 

Interestingly, we could also reverse this logic to deduce that the conservation of energy in the 
Rindler wedge forces upon us the existence of chiral edge states localized on the entangling 
surface to absorb the anomaly. In fact, we can now see that this conclusion is actually 
independent of the domain wall construction, and is true whenever a spatial entanglement 
cut passes through a topological phase, i.e., when the entanglement cut passes through a 
gapped region with a topological effective action. For example, consider the situation in 
Fig. 10 where we have two domain walls of opposite chiralities. Now let us bring the two 
domain walls closer such that they eventually annihilate — i.e., the limit where the two walls 
overlap and the low-energy modes can tunnel between each wall. In this limit there is no 
net chirality on the wall, and the domain wall states can be locally gapped. Interestingly, 
even in this situation, we are still left with chiral states along the entire entangling surface 
away from the wall, as should be clear from Fig. 10. As the wall region is shrunk to zero this 
shows that the entire entanglement cut will bind low-lying chiral entanglement modes. This 
gives us an intuitive way to understand why the low-lying spectrum of the entanglement 
Hamiltonian in a topological phase is in close correspondence to the edge spectrum in the 
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presence of a physical boundary; at the very least the low-lying spectrum of the modular 
Hamiltonian must be able to absorb the anomalous response of the domain wall/edge states. 


2.4 Explicit calculation for a Chern insulator lattice model 

We now present some explicit lattice model calculations to demonstrate the above physics in 
the context of a condensed matter system called the Chern insulator [24], The Chern insu¬ 
lator is a 2-|-l-d lattice model of massive Dirac fermions which can generate a nonzero bulk 

Chern number in certain parameter regimes. A non-vanishing Chern number Ci indicates 

2 

that the system will have a bulk quantum Hall effect with Hall conductance a^y = ^Ci, and 
thus will have chiral edge states with cl — cr = Ci. 

The 2-|-1-dimensional lattice model we consider has discrete translational symmetry with two 
degrees of freedom per site. The real-space model can be Fourier transformed to generate 
the Hamiltonian 

» = (49) 

k 

where a,b = 1, 2. The momentum-dependent matrix Hab{k) is the Bloch Hamiltonian. To 
be concrete, we consider the following single particle Bloch Hamiltonian: 

H{k) = Vx sin kxcr^ + Vy sin kya^ + {2 — m — cos kx — cos ky)a^ (50) 

where we have set the lattice constants ax, ay = 1, h = 1, kx,ky label the lattice momenta, 
m,Vx,Vy are parameters, and we assume that Vx,Vy > 0. The two bulk bands of this model 
have energies E± = sin^ kx + Vy sin^ ky + {2 — m — cos kx — cos ky^. This Hamilto¬ 

nian is a gapped topological insulator with Ci = —1 in the regime 0 < m < 2, or with 
Cl = -|-1 in the regime 2 < m < 4. Since a system with a nonzero Chern number exhibits 
localized chiral edge states, it can be used to test our hypothesis that the parity-odd change 
in the entanglement entropy of a chiral fermion is compensated by the change in the bulk 
entropy when boosted with a boost parameter r]. 

Let us now try to see how we can manipulate the parameters in this system to get something 
akin to a boost as in the relativistic cases we have been discussing in previous sections. 
Consider the system in a spatial cylinder geometry with open boundaries in x, and periodic 
boundaries in y (see Fig. 11). The dispersion of the boundary fermions localized on the x 
edge for this model is given by Vy sin ky [60]. We can thus infer that the number Vy determines 
the velocity of the low-energy chiral fermions. Indeed, the low-energy Hamiltonian of a chiral 
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fermion localized on the x-edge of the system in eq. (50) is given hy H = Vyky, when expanded 
near ky = 0. Under an inhnitesimal boost as in eq. (8), and very close to the origin in Fig. 3, 
we can calculate the change in the Hamiltonian by calculating 6^H = [H, ,^] with ^ given by 
-iv{ypt+tpy). 


kdH _dH d^\ 
\ dx°‘ dp a dpa dx°^) 
.dHdj 
dpt dt 


(51) 


where a = t,x,y. Using the fact that dH/dpt = Vy, we see that the new Hamiltonian in 
the boosted frame at hrst order in rj is given by Hboost = Vy{l + rDpy. Thus, a boost acts 
just like a change in the velocity of the low energy limit of the chiral fermions, and we will 
interchangeably use 5v/vy and p from now on in this section. 


Now let us consider the entanglement of this system. Since the bulk model is gapped (except 
at the edges) and does not have topological ground state degeneracy, it is expected to have an 
entanglement entropy which scales linearly with the length of the entanglement cut (passing 
through a gapped region), without a universal sub-leading correction: 

Sci = (52) 

where L is the length of the entanglement cut and a is a constant which will generally depend 
on the parameters in the model, and on the boundary conditions that are imposed on the 
system. 



Figure 11; We show the geometry of the system and the entanglement cut. The system is taken on a 
cylinder with open boundary conditions in the x-direction (along the axis of the cylinder) and periodic 
boundary conditions in the y direction (along the circumference). The entanglement cut is along the x 
direction and is shown in the dashed lines in the figure. The red region specifies the region on which 
the velocity parameter Vy is changed with the profile shown in Fig. 3 in the x direction, but it does not 
change as a function of y. 
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To observe the change in the entanglement entropy nnder a boost, we need to change the 
parameter Vy for the Chern insnlator (recall we assnme Vy > 0). When we have open 
bonndary conditions, and with the cnt nsed in Fig. 2, we expect there to be no anomalous 
change in the parity-odd contribution to the entanglement entropy of the system under 
a uniform change in Vy (i.e., a, boost). In this case, since both edges will be boosted 
identically, then the sum total of the change in the anomalous parity-odd contribution to 
the entanglement entropy should vanish. We do note, however, that the total entanglement 
entropy may still be modified when changing Vy due to, for example, changes of the non- 
universal parameter a in eq. (52). 

Instead, to isolate a single edge as in the domain-wall configuration, we change the velocity 
parameter Vy near only one of the edges of the Chern insulator lattice system, and have it 
return to its initial value deep in the bulk, but far from the other edge^^. In this case, if there 
is to be no anomalous entanglement contribution, then the change due to the single edge must 
be compensated by the boosted bulk. An important difference between the calculations we 
do in this subsection versus the rest of the paper is that we perform global boosts on the edge 
(as done in Ref. [28]). We focus on such global boosts here to maintain translation invariance 
in the lattice along the edge direction which simplihes the calculations considerably. The 
expected answer for the change in the edge-state entanglement entropy under a global boost 
is for a single interval with two endpoints on the edge. This has an extra factor of 

two compared to the answer that we found under the local boost prescribed in Fig. 3 and 
eq. (44). 


Before we calculate the bulk-boundary cancellation, let us hrst show that a lattice chiral 
fermion itself will have a non-vanishing change in the entanglement entropy when boosted. 
For free-fermions the correlation matrix method [61] is a convenient technique that we will 
use to calculate the change in entanglement entropy of a chiral fermion under such a boost. 
Now, we could take the full bulk spectrum of our 2-|-l-d lattice model, project onto the 
chiral edge modes, and then calculate the entanglement of this subset of states. However, 
for clarity of presentation, and calculational simplicity, we note that the chiral edge modes 
are essentially identical to the subset of right-movers (or left-movers) of a l-|-l-d single-band 
tight-binding model of lattice fermions. A hopping model which realizes this analog is 


H 



1 


c]c,+i 


(53) 


which has a dispersion E{ky) = VySinky. This model mimics the lattice chiral fermions 

is sufficient to change it on a region extending into the bulk from the edge that is much larger than 
the penetration depth of the localized edge states. 
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except that there are both right and left movers since the dispersion is dehned over the 
entire Brillonin zone. This is a conseqnence of the Nielsen-Ninomiya theorem. Fortnnately, 
since the model is non-interacting, we can isolate jnst the left or right movers and calculate 
their entanglement using this simpler model. 

To calculate the entanglement spectrum and entropy we need to construct the correlation 
matrix. The correlation matrix for this hopping model in l-|-ld is given by [61] 

(54) 

where n{qy) is the occupation number of the “edge states” as a function of qy, and we 
have used the matrix element (cj^Cj) = . The calculation and end result, for the 

entanglement entropy is discussed in detail in Appendix C. We hnd that varying the entropy 
with respect to the cutoff momentum /cq is ^ x tan Here /cq is the momentum 

cut-off for the occupied states in the chiral fermion branch, and it is of interest in this 
calculation, as shown in Appendix C, since changing the velocity of the chiral fermions acts 
identically to changing the momentum cut-off. Interestingly, we hnd that our result for the 
boosted entanglement entropy has a lattice factor tan^, which is precisely related to the 
lattice factor appearing in the Calabrese-Cardy lattice formula for the spatial entanglement 
entropy of a CFT [1]. At a cutoff momentum of fco = f; we recover the exact continuum 
value for a global boost on the edge. When ~ 0,7r, we expect the lattice effects to be 
maximal because we are close to the limit of an empty or fully hlled edge band respectively, 
and the dispersion becomes non-chiral. The fact that we have found lattice modihcations 
to the connection between entanglement and gravity may give hints about how to properly 
discuss lattice regularizations of gravitational anomalies, although we leave that to future 
work. 

We have now conhrmed that a free chiral fermion gives a modihcation to the entanglement 
entropy when boosted. When the entire bulk and boundary system is considered together 
(with open boundary conditions in x and periodic boundary conditions in y as illustrated in 
Fig. 11), we expect the anomalous change in the parity-odd contribution to the entanglement 
entropy to cancel out. To extract the parity odd part of the entanglement entropy, our 
strategy will be to consider the system at positive and negative (which correspond to 
positive and negative Chern number respectively, i.e., opposite parity phases) and change 
the velocity by a small amount v± = —>■ + 6v near one of the edges to isolate the 

change from the other edge. The parity odd part of the change in entanglement entropy is 
then calculated as 

t = • ( 55 ) 
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Figure 12: We plot the change in entanglement entropy 6S vs the length of the entanglement cut 
(the rectangle including contributions to the length along each edge, and from the sides of the cylinder) 
of a Chern insulator system with open boundary conditions in the x-direction and periodic boundary 
conditions along the y-direction with the entanglement cut from Fig. 11 and the boost near a left moving 
edge. The blue cross markers are the simulation results and the yellow line is a linear fit. We notice 
from the inset fitting equation that the plot is linear in the length of the entanglement cut and the tiny 
constant offset is likely from finite size errors. The anomalous contribution from a chiral edge state 
would show up as a finite offset on the linear fit of order = —g), which we do not see. This 

indicates cancellation between the bulk and boundary contributions to the entropy. The parameters 
used in the plot were Vy = 1.0, 6v = 0.01 and m = 1.0 from the Chern insulator model in eq. (50). 

We have seen that a change in entropy due to the boundary chiral fermion should yield 
a constant contribution to (or equivalently ^) which is independent of the size of the 
entanglement cut. However, what we find in Fig. 12 is that the total change 6S is perfectly 
linear in the length of the entanglement cut, and does not have any signature of a constant 
offset coming from the chiral edge state. In fact, we can understand the linear-scaling piece 
as follows. The entanglement entropy from the bulk Sbuik{vy^^) = a{v^^)L is an even function 
of the velocity Thus, in the small 6v regime, there is an additional non-universal change 
in eq. (55) coming from a{v^^) that can be approximated as ~ a\vl^^)L + 0{6v). This 
is where we expect the non-vanishing contribution in Fig. 12 comes from. 

From the linear fitting equation shown in Fig. 12 we see that the constant offset is essentially 
vanishing, and hence we have demonstrated that the bulk variation does indeed cancel off the 
universal part of the variation of the boundary entanglement entropy as expected from our 
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earlier calculations. Additional consistency checks were done to check that the plotted in 
Fig. 12 is independent of the rj = 6v that was chosen in the small r]/v^^ regime. The result 
was also found not to depend qualitatively on the the parameter m from the bulk Chern 
insulator model in eq. (50), as long as one remains in a gapped topological phase. 


3 Discussion 


In 2+1 dimensional parity violating theories, domain walls between topologically distinct 
phases support localized 1+1 dimensional gapless excitations. The low energy effective 
description of these excitations is in terms of (generically) chiral conformal held theories, 
which possess a variety of anomalous currents. In this paper, we considered the entangle¬ 
ment entropy in such 2+1 dimensional theories, in the situation where the entangling surface 
intersects the domain wall. We have shown that the anomalous dependence of the entangle¬ 
ment entropy on local boosts of the entangling surface in the 1+1 dimensional theory are 
accounted for by inhow from the bulk. We have thus termed this phenomenon entanglement 
inflow. This is in keeping with the general structure of anomaly inflow, and is a manifestation 
of the fact that the 1+1 dimensional theory does not really exist in isolation. It is satisfying 
that entanglement entropy can be consistently dehned (i.e., independently of local boosts 
which preserve the entangling surface) by going to the full bulk description. In a manner of 
speaking, the bulk theory provides a UV completion of the edge CFT, and it is only natural 
that one has to go to the bulk to dehne entanglement entropy in a consistent way. 

Interestingly, we were able to give a physical picture for cancellation of the bulk and bound¬ 
ary entanglement responses in terms of (entanglement) energy inflow in the Rindler wedge, 
mediated by the low-lying chiral states in the entanglement spectrum. We inferred the 
existence of these states from the requirement of anomaly cancellation, or in other words 
(background) diffeomorphism invariance, in the Rindler wedge. Turning the logic around, 
one can deduce the existence of low-lying chiral states in the entanglement spectrum from 
the physical principle of energy conservation in the Rindler wedge, through thought experi¬ 
ments like the one considered in Fig. 10. This gives an interesting new vantage-point to think 
about the correspondence between the low-lying spectrum of the entanglement Hamiltonian 
and the spectrum of physical edge states in topological phases. 

We end with some discussion of open problems and future directions. First, it is interesting 
to consider similar phenomena in higher dimensions [29]. Gravitational anomalies exist in 
4/c + 2 dimensions (for k = 1,2---), and we expect our arguments to carry through in 
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these cases as well. Of course, in order to study these higher dimensional analogues, one 
presumably has to admit more complicated geometries for the entangling surface. Second, the 
effect of parity-violation on entanglement entropy is by itself an interesting arena to explore. 
Another interesting question is whether the above physics can be described holographically in 
strongly coupled held theories. Indeed, the anomalous response of the entanglement entropy 
to boosts was considered in the context of holographic entanglement entropy in [28], and more 
recently in [62,63]. In our case, in order to describe a holographic CFT with three spacetime 
dimensions and a parity domain wall which hosts a lower dimensional chiral CFT, one might 
consider asymptotically-AdS Janus-like solutions [64] supported by pseudoscalars [65]. It 
might be an interesting exercise to consider the Ryu-Takayanagi entanglement entropy in 
this context [66] . Since most of the arguments in this paper were based on general principles, 
there is no reason to suspect any radical modihcations at strong coupling. Additionally, 
it would be interesting to consider the properties of other symmetry-protected topological 
phases which are not chiral, but still harbor low-energy edge states and entanglement modes. 
We leave these questions for future exploration. 

Note: While this paper was in its hnal stages of preparation, references [30,32] appeared 
on the arXiv, which have some overlap with the contents (in particular, section 2.2) of this 
paper. 
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Appendix 


A Perturbative change in Entanglement Entropy 


In this appendix, we want to prove the formula (11) for the perturbative change in en¬ 
tanglement entropy. Consider a generic quantum field theory that admits a path-integral 
description in terms of a field 0, which collectively denotes all the fields over which we will 
integrate. The density matrix corresponding to the ground state wavefunction at f = 0 is 
given by 

e-0 

p = lim/3^oo -7-TV (56) 

Tr^ 

where H is the Hamiltonian, and Ti is the relevant Hilbert space. In the path integral 
language, we can describe a matrix element of this density matrix as a product of path 
integrals over the regions te > 0 and < 0 of Euclidean space with the appropriate 
boundary conditions 

(0-|p|0+) = 4/ [D(p]rE<oe~^^^^ [ [T>0]^^>o (57) 

^ ./(/)(0-,x)=</._(x) d0(O+,x)=(/>+(x) 

where we have denoted the spatial coordinates collectively as x = {y,x), and Z is the 
partition function of the theory on 


Let us now denote the reduced density matrix for the half space A = {x^ = {y,x)\y > 0} by 
Pa- The matrix element (0 ^)|pa| 0+) is then given by gluing the above path integrals along 
the complementary space H at rg = 0 

(</>(!Ipa 10+) = 7 / [^0] (58) 

where x"^ are spatial coordinates on A, and 0+ denotes a field configuration on A. By slicing 
this path integral along the angular direction p in the {te, y) plane, it becomes immediately 
clear that the reduced density matrix can be written in operator form as [67] 


PA 


e-2nK 


(59) 


where K is the generator of (p-rotations 


K 



dyyT {0,y,x). 


(60) 
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Equation (59) shows that, up to an overall shift coming from the normalization, the entan¬ 
glement Hamiltonian in this case is given by 

He = -In PA = 2ttK + In Tr^^ . (61) 


Now, if we turn on a small (background) metric deformation we hnd 
( 0^1 {pA + Spa) 10 +) = 7 -^ / [E>0] + ^ drEd<^-^x 5gf^-iTE,x)T^ArE,x)+- 

[Z + dZ) y^(o+,x^)=</.^ 

(62) 

and therefore to linear order in 6g, we have 

SpA=^f [D(j)]e~^^S] f dTEd‘^~^yi6g^‘'{TE,:>i.) (%u{rE,^) - ■ (63) 

2Z J^R Jr3 V / 


From equation (59), we can then infer the following operator expression for the change in 
the reduced density matrix 

SpA = 7^ [ 6g^''{TE, x) ('pA^0j.(r, x) - pA{f^y{T, x))) (64) 

The operator T{te,^) above is to be interpreted (from the point of view of the reduced 
density matrix) as a Heisenberg operator 

f>^^{p,r,x) = {R-\^))\{R-\^))\e^^f^-{0,r,x)e-^^ (65) 

where {p, r) are polar coordinates in the (rg, y) plane, and R{p) is the appropriate rotation 
matrix. 

Finally, from the dehnition of entanglement entropy, we have 

SSee = Tr^^ {SpA He) (66) 

Putting everything together, we hnd 

SSee = t; f dTEd'^~^^Sg^''{TE,^)Z^T^HA(pAT^,u{T,^)HE-P a{T^^{t,^))Hj^ 

^ jR3 V / 

= ^ j ^dTEd'^~^yi5g^''{TE,^)(^(T^^{T,yi)HE) - (He)) 

= aI dTEd'^~^yi6g^''{TE,yd) (T^y{T,yi)HE) (67) 

Z J^Z \ / conn. 


where the correlation functions appearing above are Euclidean correlation functions on M^, 
and in the last line we have used the dehnition of a connected correlation function. The 



above integral is divergent in general, and one must be careful to regulate it by removing a 
tubular neighbourhood around the point te = y = ^ 



dTEd"^ 6g>^''{TE, x) x)He'^ 


conn. 


( 68 ) 


Thus, we arrive at the expression used in the main text. 

B Entanglement spectrum for free Dirac fermions 

In this section, we demonstrate explicitly that in the case of a free d = 2 + 1 massive 
Dirac fermion with a mass domain wall, there are chiral modes along the domain wall, and 
more importantly, along the entangling surface. We will work in Lorentzian signature in 
this section, because the spectrum of the entanglement Hamiltonian will be much simpler 
to interpret in this case. In order to get the spectrum, we simply write down the Dirac 
equation on the Rindler wedge in with a domain wall prohle for the mass m{x). We 
then conformally map to M x (this is not necessary, but we do so to make contact with 
the discussion in the main text), and solve the resulting Dirac equation on M x 

The Dirac equation on a general background is given by 



(69) 


where is an orthonormal frame on spacetime, u°‘b is the spin-connection, and 7 “ are 


matrices satisfying the Clifford algebra { 7 “, 7 ^} = 27 “^!. A simple property of the Dirac 
equation is that it is covariant under the transformation [ 68 ] 


e“ = De“, T = D 2 T, m = D 


(70) 


where D is an arbitrary function on spacetime and d is the spacetime dimension. We can use 


this fact together with the conformal equivalence of the Rindler wedge with x to solve 
the Dirac equation on the latter geometry with the mass m = Q~^m. Let us parametrize 
M X with the coordinates (T, z, x) with T being the time. It is convenient to use the 
following coframe 



The (Levi-Civita) spin connection in this case is given by 



( 72 ) 
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and the mass is given by m = zm{x). The Dirac eqnation on x is (dropping the tildes 
for convenience) 

- |) + 'y'^zdx + zm{x)^'^ = 0. (73) 

We’re interested in the sitnation where m{x) has a domain wall prohle, which means that 
m{x) interpolates between —mo and +mo as x runs from —oo to +oo, mo being a positive 
constant. For simplicity, we will take m{x) = tuq sign(a;). 


In the presence of the above mass domain wall prohle, one hnds a special normalizable 
solution to the Dirac equation of the form 

Tg;^(T,x,a;) = (74) 

uj{k) = k, 7Vo(/^) = +i’o{k) (75) 

This is a zero-mode of the Dirac operator which is localized on the domain wall between the 
two-phases m < 0 and m > 0 (hence the subscript DW), and is chiral with respect to the 
domain wall chirality matrix 7 ^. In the limit m —)■ cxd (i.e., at energy scales small compared 
to the mass), we can model the dynamics of this mode as a 1 -|- 1 -dimensional chiral fermion. 


For |a;| >> 0 (i.e., far away from the domain wall), it is a simpler matter to construct generic 
solutions. Going to Fourier space in T and x 


^{T,z,x) = r ^ r 


27r 


— CKD ^ —OO 


27r 


(76) 


and writing the Dirac equation in terms of '0(±) = |(1 ± 7^)4/, we obtain 

izdzi>(+){uj, k] z) - {u + 1 ) '0(+)(a;, k; z) = z{ik - t, z) (77) 

izdzil)(_){ijj, k] z) + [u - 1 ) k] z) = z{ik + m)'0(+)(a;, k] z) (78) 

Eliminating '!/’(-) j dehning '0(+) = X'0(+) we obtain 

z‘^dl'4)(+){u:,k-,z) + zdz'ip(^+){uj,k-,z) - [\ + iu)^ + z^+ k‘^) V^(+)(a;, fc; z) = 0. (79) 

Therefore, generic bulk solutions can be written in terms of modihed Bessel functions^^ 
V’(+f ^(T, z,x) = z Ki^-^ (VF + m2zj ^^uT+^kx 


l(UI,K)/rp \ k -\- iTfl 

^(_) {T,z,x) = -j===z 


Vk‘^ + 




m2^) Xo{uj,k) 


( 81 ) 
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We have dropped solutions involving + m?z) because they diverge for large z. 
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which are oscillatory for small z and exponentially decay for large z. To identify modes 
localized on the entangling snrface, care mnst be taken to regnlate the theory. Here, we 
choose to cnt off the Rindler wedge at z = a, with “brick-wall” bonndary conditions [50] at 
z = a 

7 ^T(T, a, x) = ctT(T, a, x) (82) 

where a = ±1. We will only constrnct solntions for the region z > a, and the above bonndary 
condition ensnres that there is no flow of charge or energy across the cntoff snrface. The 
bonndary condition implies 

Ki^i^ (aVk'^+ 771"^^ = (a^FT^) . (83) 

This eqnation can be solved nnmerically to obtain a; as a fnnction of k, namely the entan- 


OJ 



Fignre 13: The entanglement spectrum corresponding to the boundary condition with a = sign(m) 
for m = 10, and a = 0.01. 

glement spectrnm. For a = sign(m), we find spectrnm (see Fig. 13) 

u{k) = vk + ■ ■ ■ (84) 

where the ellipses denote terms of order (^)^. The velocity v in this case is non-nniversal, 
and is given by 

p — 2ma 

V = -;- r. (85) 

2mF(0, 2ma) 

For \k\ small compared to \m\, the wavefnnctions are localized near the regnlarized entangling 
snrface z = a, and are chiral with respect to the entangling snrface’s chirality matrix 7 ^, np 
to 0{k/m) corrections. Thns, we see that the cutoff entangling surface hosts a chiral mode. 
Quantizing this mode gives us Casimir momentum density on the entangling surface, which 
is what was used in the main text. 
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The other branch a = —sign(m) gives rise to gapped states in the entanglement spectrnm. 
These states are oscillatory for 2 ; large compared to a, and shonld be interpreted as bnlk 
states. Since the entanglement entropy is dominated by the low-lying spectrnm of the entan¬ 
glement Hamiltonian, it snffices to focns on the low (entanglement) energy states, which in 
the present case consist of the chiral modes localized on the domain wall and the regnlated 
entangling snrface. 


C Correlation matrix and entanglement entropy for an 
edge chiral fermion 


The correlation matrix as given in eq. (54) is 





dqy 

27r 




( 86 ) 


where n{q) is the occnpation nnmber as a fnnction of momentnm. For the case of a right 
moving edge state with the dispersion E = Vy sin qy which, for a single chiral fermion does 
not traverse the entirety of the edge BZ, n{qy) = 1 for —ko < qy < 0 (assnming > 0) 
where /cq is a cntoff determined from where the edge band Vy sin qy meets the bnlk band E_. 
We can thns simplify the correlation matrix to 




1 - e ^ ^-iko(m-n)/2 ^^^ko{m - n)/2 

27ri{i — j) 7r{m — n) 


(87) 


The eigenvalnes of this matrix when i, j are restricted to a snbsystem can be nsed to calcnlate 
the entanglement entropy. To move from a right moving edge to a left moving edge, we simply 
need to flip ko —)■ —ko in eq. (87). This corresponds to choosing n{qy) = 1 for 0 < qy < ko- 


The phase factor from (87) does not affect the eigenvalnes of the correlation matrix, which 
we can see as follows. The eigenvalne eqnation for the correlation matrix is given by 
Cij(f>k{j) = Ck4’k{i)- Hence, the phase factor in Cij can be absorbed into the 4’kUy^ 
so that C,k remain the same. Modnlo the phase factor, the correlation matrix from this prob¬ 
lem is identical to what is fonnd for a ID tight-binding electron, bnt with a Fermi wavevector 
of ko/2 [69]. The entanglement entropy is given by 




log(l + e "*)-!- 


efc 


1 -f e’^'= 


( 88 ) 


where are the eigenvalnes of the entanglement Hamiltonian and are related to (k as 


(k — 


1 


e‘k + 1 ■ 


Qnoting from Refs. [61,69], we have that ^ = 


log 2L sin ko/2 


where L is the 
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length of the subsystem in question. This can be used to perform the integration in eq. (88), 

S = ‘^^°s‘iLsmko /2 j^oo |^iog(l + 6“'^'=) + where each term in the integration gives us a 

2 

factor of There is an extra factor of two overall which we have inserted when converting 
from the sum to the integral to account for both positive and negative e^. The end result is 
that S = ■ Under a variation of the cutoff we see that 


5S ko 1 

—— = tan — X - 
6ko 2 6 


(89) 


for a right moving edge. So, in general we have ^ ~ tan ^ x 
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